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Abstract 



We discuss spherically symmetric, static solutions to the SU(2) a-model on a de 
Sitter background. Despite of its simplicity this model reflects many of the features 
exhibited by systems of non-linear matter coupled to gravity e.g. there exists a 
countable set of regular solutions with finite energy; all of the solutions show linear 
instability with the number of unstable modes increasing with energy. 



1 Introduction 



With the unexpected discovery of spherically symmetric static solutions to the Ein- 
stein Yang Mills system by Bartnik and McKinnon [||] new interest focused on 
particle-like solutions for non-linear matter coupled to gravity. Many of these sys- 
tems under consideration have been shown to allow for asymptotically flat soliton 
and black hole solutions (for a recent review on the EYM-system see [0). More 
recently Brodbeck et. al. || investigated the EYM-equations including a cosmolog- 
ical constant and showed, that a rich spectrum of solutions exists depending on the 
dimensionless parameter in the theory. 

An extremely simple model that lives on a compact 3-space was previously consid- 
ered by Bizon [Q: he looked at harmonic maps from (S 3 ,G) to (S 3 ,G), where G 
denotes the metric of constant curvature on S 3 . In the present paper we modify his 
model to harmonic maps from de Sitter spacetime to S 3 , i.e. we consider a SU(2) a- 
model that lives on a de Sitter spacetime with positive cosmological constant. This 
model has the advantage of being more physical while keeping the simplicity. We 
look for spherically symmetric static solutions in a static coordinate patch of this 
geometry. Although the slices of constant time are S 3 , the field equations resulting 
from this model differ from those considered by Bizon. Numerical analysis of these 
equations reveals the existence of a discrete sequence of static solutions with finite 
energy. The energy of this sequence is bounded from above by a limiting (singu- 
lar) solution. We show that all of these solutions are unstable under linear (time 
dependent) perturbations, the number of unstable modes increasing with energy. 
Some care is taken to analyze the high excitations. To check the eigenvalues of the 
unstable modes we compare them to the perturbations of the limiting solution. For 
this we use a self-adjoint extension of the operator governing these perturbations. 



2 The Non-linear cr-Model 

We take (M,g) to be de Sitter space (with positive cosmological constant A). The 
non-linear cx-model under consideration is a harmonic map X : (M,g) — > (S 3 ,G), 
defined by the stationary points of the action Sm = / d A xC,M with the Lagrange 
density 

Cm = -V^g^gTV , ii X a V v X b Gab{X). (1) 
This yields the field equations 

(V M V,X A + t A c V^X B V v X c ) = 0, (2) 

where T^ c denote the Christoffel symbols with respect to the metric Gab of the 
target manifold. We choose coordinates (X A ) = (/, 6,$), so that G assumes the 
form Gab = diag(l, sin 2 /, sin 2 / sin 2 6). It is well known that part of the de Sitter 
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space can be covered by a static frame i.e. coordinates (t, r, tp) in which the metric 
becomes 

ds 2 = -(l-^\ dt 2 + (l - ^ 1 dr 2 + r 2 d$ 2 + r 2 sin 2 $dp 2 , (3) 

with < r < a and a 2 = 3/A, A being the cosmological constant. The singularity at 
r = a signals the breakdown of the static coordinates. We seek for regular solutions 
of Eqs. (0) in the domain, where the metric takes the static form (|3]) and require 
the mapping to be static, i.e. dtX A = 0, and spherically symmetric. The latter is 
guaranteed by the Hedgehog- Ansatz: 

f(af) = f(r),e(a?) = ( 4 ) 

In order to map the pole r = to a single point of S 3 , we require the function f(r) 
to obey the boundary condition /(0) = 0. Inserting this Ansatz into Eqs.(@) yields 
the non-linear ordinary differential equation 

r2 ( 1 "5) / ') - sin2 / = - ( 5 ) 
The energy associated with an excitation of the form (|4|) is given by 

E = 4tt/ 2 f% 2 dr{{\ - 4)f 2 + ^sm 2 f} (6) 

where a factor of 2 has been included, because the r-integration runs twice from 
to a. 

Note, that the energy functional E[f] ([]) may be interpreted again as a spherically 
symmetric harmonic map, this time from (S 3 , (1 - 4) 2 G) to (S 3 , G). 

To our knowledge the only analytic solutions to @ are the trivial solution / = 
(mod 7r) with energy Eq = and the solution / = n/2, which will be called singular, 
since it fails to obey the boundary condition at r = (see section H). 



3 Numerical Results 

In the following we introduce a new independent variable x defined by r = a sin \ 
and x = — lntanx/2. The coordinate x tends from x = (r = a) to x = oo (r = 0). 
We also define a new dependent variable h by h = f — tc/2. In these new variables 
equation @ becomes 

h" - ftanh x — — —) h' + sin 2h = 0, (7) 

V sinh x cosh x J 

where prime now denotes differentiation with respect to x. 
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Equation (|7|) is similar to the one studied by P.Bizon [[|. Ours differs from his by 
the term (sinhxcosha;) -1 which causes a singular point of the differential equation 
at x = (corresponding to the breakdown of the static frame at r = a). For the 
solutions to be regular at x — 0, h'(0) has to vanish. The other singular point is at 
x = oo where h(oo) = ±tt/2 (mod7r). The formal power series expansions of h(x) 
around these singular points give: 

h(x) = b — l/4(sin26)x 2 + 0(x 4 ) for x ~ and , , 

±h\x) = n/2-2ce~ x + 0(e~ 3x ) for x -> oo. {b) 

One can show that near x = as well as x = oo, there exists a family of solutions, 
each of them depending only on a single parameter, b = h(0) and c respectively. 
Since with h(x) also h(x) = nix ± /i(:r) is solution to equation (|7|), we restrict the 
parameter b to the range < b < ir/2. 

We integrated equation (|7|) numerically, starting at x — 0, trying to find values for b, 
such that h(x) shows the correct asymptotic behavior for large x. With the values 
found for b we used the two-point shooting and matching method to determine 
the parameter c. Both numerical routines, d02cbf and d02agf, were provided by 
the NAG-Library [H]. The numerical analysis suggests the existence of a countable 
family of globally regular solutions of Eq. (|7]) which may be labeled by their number 
of zeros, and indeed a rigorous proof of existence can be given along the lines in • 

Figure p] shows the first four solutions h n (x), where the sign of b n has been chosen 
such that h n {oo) = —n/2. The numerical results b n and c n for the first four solutions 
are displayed in Table |]. 



n 


b n 




b 2 lb 2 
u n/ u n+l 






1 


0.4458536 


5.0270678 


10.636094 


10.831606 


1.8304231 


2 


0.1367103 


54.451223 


10.734453 


10.747914 


1.9842302 


3 


0.0417264 


585.23706 


10.747677 


10.748889 


1.9985328 


4 


0.0127278 


6290.6484 


10.748955 


10.749068 


1.9998635 


'oo' 






10.749087 


10.749087 


2 



Table 1: Numerical values of the parameters b n , c n , the ratios 6 n /6 n+1 ,c n+ i/c n and 
the energy E n of the first four excitations. 



4 Large Excitations 

In what follows we give arguments for understanding the numerical value of the 
ratios formed from the integration constants for consecutive values of n as listed in 
Table [l]. The line of reasoning, which is based on approximating Eq. (|7|) for large 
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Figure 1: The first four solutions h n (x) (The sign of b n has been chosen such that 



ft n (oo) 



excitations by linear respectively translation invariant differential equations, is the 
same as in ||, but for completeness we will repeat it here. 

For large n one can distinguish essentially three regions for equation (0): 

Region I, for < x < 4, where numerical analysis tells us, that ft is small, so that 
Eq. (0) is well approximated by 



V smh x cosh x 



ti + 2ft = 0. 



(9) 



Since equation @ is linear ft„(a;) ~ b n g(x) for large n, where solves equation 
(H) with initial conditions p(0) = 1 and </(0) = 0. 

In region II, 4 < x < 2n, where ft. is still small, one can approximate equation (|^) 
for large n by 

h" - ti + 2ft = 0, (10) 

which is linear as well as translation invariant. Eq. (1Q|) is solved by ft n ~ b n g(x) 
with ^(x) oc e^ 2 sin {y/lx/2 + <5). 

Finally in region III, for x > 2n, equation (|^) becomes 



ft" - ft' + sin 2ft = 0, 



fill 



and is solved asymptotically by ±ft n (x) = 7r/2 — 2c n e x . 
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Now, by using translation invariance of Eqs. (|I0D and ( |TT| ) as well as the asymptotic 
behavior of h n , one obtains in region II and III for two consecutive solutions: 

h n+ i{x) ~ h n (x - lnc n+1 /c„), (12) 

if we choose h n+ i and h n to have the same asymptotic behavior, (which is possible 
by choice of sign). In the regions I and II, due to the linearity of Eqs. ([D and fllO|) 
and h n (0) = 0, h n+ i and h n are connected by 

K+i = -r^-h n . (13) 

On 

Using the solution to Eq. ( |i~0|) and the fact, that, in region II, h n+ \ has one zero more 
than h n , one determines the parameter values of consecutive solutions to behave like 

bl/b 2 n+1 ~ c n+1 /c n ~ e 2 ^ (14) 

for large n. 



5 Stability Analysis 

We now discuss linear stability of the solutions h n under time dependent spherically 
symmetric perturbations of U, i.e. we use the Ansatz: 

h(r,t) = h n {r) + Sh(r,t) (15) 

where h n (r) is solution to the static field equation ([?[) and Sh(r,t) is considered to 
be small, so that we retain only terms linear in Sh. The full time dependent field 
equation is: 

( r2 i 1 ~ 5) h ) + sin 2h = (1^4) k (16) 

Inserting fllS]) into ( p^) and linearizing in Sh gives: 



2 \ \ ' 2 



r 2 1 - — \5h'\ +2 cos(2/i n )5/i = ^5h. (17) 



a 2 



We set 

5/i(r,t) =£(r)e- iff *, 

which leaves us with 



r 2 (l-^Vl +2cos(2^)e= t: ^ 2 \ e (19) 



a 2 
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If this eigenvalue problem has any modes with a 2 < , 5h would grow exponentially 
in time and our solutions would be linearly unstable. 

The energy of the time-dependent field is given by 

E = 4./| j° r\lr [jp^f + (' ~ 9} h " + ? ^ *) 

Regular perturbations have to satisfy the boundary conditions £(a) = and £(0) = 
0. Using the coordinate x Eq. fll9"|) reads 

sinh 2 (x)£" + tanh(x)(l - sinh 2 (x))^' + 2cos(2/i n ) sinh 2 (x)£ = -A£, (20) 

where A = <t 2 q; 2 . 

The asymptotic behavior of £ is: 

£ ~ ,x v/ ~^ for x ~ and 
£ ~ e _x for x — > oo. 

Using routine d02kdf of the Fortran NAG-Library we find numerically (up to 
n = 3) that the n-th "excitation" has n negative eigenvalues, i.e. n unstable modes. 
The numerical values of the first three negative eigenvalues of Eq. fl20D for the first 
six excitations are displayed in Table [2| 



n 


A? 


K 


A3 


1 


-2.111967 






2 


-2.065756 


-183.16864 




3 


-2.061672 


-180.59393 


-21007.87 


4 


-2.061294 


-180.3407 


-20715.26 


5 


-2.061270 


-180.3169 


-20686.48 


6 


-2.061256 


-180.3148 


-20683.78 



Table 2: Numerical values of the first three negative eigenvalues 



6 The Limiting Solution 

From Eq.([L4]) we conclude that the amplitudes in region I and II decrease with n 
as e~ n7T ^ . At the same time the region of oscillations extends to larger and larger 
values of x. Thus for n — > 00 the solution h n (x) tends point wise to zero for any 
finite x. 

The limiting solution h^ix) = does not satisfy the required boundary conditions 
at x = 00, and as a consequence its energy density is singular. Nevertheless the 
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total energy (||) is finite, Eoo = 8nf 2 a, which can be shown to be the upper bound 
for the energy of all regular solutions h n : 

Expressing the static Eq. (||]) in terms of h, multiplying by h leads after a partial 
integration to the relation 



„2 



drr 2 (l -)h'= drhsm2h. (22) 



„2/< '",,,-> 


Thus for regular solutions the energy functional E[h] can be written as 

E[h] =4vr/2 J dr(h sin 2h + 2 cos 2 h). (23) 
o 

Since the regular solutions satisfy \h\ < it/2, the integrand of ( p3|) has a maximum 
for h = 0. 

The reason for our interest in the singular solution h(x) = is, that the eigenvalues 
of the linear perturbations about regular static solutions seem to converge to some 
value A£° for fixed k as n — > oo. (See Table ^). So the question is, can we find the 
limiting values by calculating perturbations of the singular solution? This would 
yield a check of the numerics and moreover tell us, that we really have obtained all 
eigenvalues in the stability analysis. 

Let us first have a look at Eq. fl20D in the limit n — > oo, i.e. = 0. Introducing 

a new variable w = — sinh 2 a; and setting £(w) = {—w) y ^'^ 2 y{w) leaves us with a 
hyper-geometric differential equation 



w{l-w)y"{w) + (^\+l- w (^/\X\ + l)jy'(w) + 



+ | J ^-^ 1 -^H = 0, (24) 

which is solved by the hyper-geometric Function F(a, b, c; w) = S^Lo ^f^"" 07 ™' 
where (a) n = a(a + 1) - • • (a + n — 1); (a) = 1 etc. and 

a = -1/4 + ^/2 + ^77/4, 

b = -1/4 + y/\X\/2 - iV7/4, 

c = 1 + y^Aj. 
Analytic extension of F(a, 6, c; w) for |iu| > 1 gives 

£(w) = (-w) 1/ ^ i ^ /i P(X)F(a,l-c + a,l-b + a;l/w) + 

+ (-w) 1/i+iV1/A P{\)F(b, 1 - c + 6, 1 - a + 6; (25) 
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where denotes complex conjugation and 



P(X) 



r(c)r(6 - o) 
r(6)r(c - a) ' 



In order to get a spectrum of eigenvalues, that may serve as a limit to the spectra of 
regular perturbations, we need to specify boundary conditions at w = — oo (r = 0). 
This is not straightforward, since the solution hoc = is already singular. In what 
follows we make use of an idea which is due to P. Bizon H: 



We first note, that the differential operator appearing in equation (2D), rewritten in 



the coordinate w 



sinh x, 



4iu(l - w) 3/2 



w 



d 



dw 



w) 1 / 2 dw / 



2w cos(2/i r 



is formally self-adjoint with respect to the scalar product 

ro dw 



i.e. 



-w)(l — w) 3 / 



A(-w) 



w 



,1/2 



(26) 



(27) 



:28) 



The boundary terms vanish for (regular) perturbations of regular solutions h n . It 
is this property that we would like to impose on the singular perturbations. For 
solutions of the eigenvalue problem ( f[lf) this requirement yields the condition, that 
the ratio of coefficients in (|25| ) p = e~ 2iarg( - p ) has to be independent of A. In fact, 
admitting solutions that behave for w — > — oo as in fl25|) with p fixed, corresponds 
to a self-adjoint extension of the operator Too, (see Appendix). 

We are now free to choose one of the eigenvalues, say = —2.06126. The spectrum 
of ^ then is determined by e 2iarg ( p ^ Afc ^ = e 2 * arg ( p ^ Al ^, or equivalently 



arg(P(A fc )) = arg(P(A0) + (k — l)vr. 



(29) 



We obtain e.g. A^ 



-180.315 and A|f 



-20683 in good agreement with Table 0. 



7 Outlook 

In this paper we have analyzed spherically symmetric solutions to the non linear 
sigma model on de Sitter space with a positive cosmo logical constant. These solu- 
tions are regular everywhere in a static coordinate patch of this spacetime i.e. up to 
the cosmological horizon. Despite the simplicity of the model it shares many of the 
properties of matter models coupled to gravity and may therefore help to understand 
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why such excitations exist. Furthermore we expect that when the back reaction is 
taken into account, at least for weak coupling, these solutions will persist. 

The authors would like to thank P. Bizon, S. Husa and H. Narnhofer for helpful dis- 
cussions, moreover we wish to thank S. Husa for his support regarding the numerical 
computations. 



A Appendix 

The perturbation equations (IT5) can be transformed into a one- dimensional Schrodinger 
equation by introducing new dependent and independent variables: 

£ = — and p = - In — ^— (30) 
r 2 a — r 

Eqs. ( |T9"D then become 

d 2 (cos(2fr n ) + tanh 2 p) 

— ~r^ ( t ) ~ 2 —2 = A0. (31) 

dpr sinn p 

For small p the potential behaves as \ for regular perturbations and as — \ for the 
singular case = 0. 

In order to determine domains (if possible) on which the operators are self-adjoint, 
it is sufficient to consider the operators r = — + \ for the regular case, r 0iOO = 

— cT* ~ T f ^ or ^ ne singular case respectively. These operators have been studied by 
H. Narnhofer 0: 

It turns out that the operator tq is self-adjoint on the domain 

V{t q ) = {0, 0'abs. cont., 0" G H, 0(0) = 0'(O) = 0}, (32) 
where H = L 2 (dp, (0,oo)). 

To.oo is hermitian but not self-adjoint on the above domain, since dim(Ker(rQ OC =p 
i)) = n± = 1. Solutions of (r^^ =f i)<fi± = are 0+ oc p 1/2 (^ iv ^/ 2 (p ei7r/4 ) - 
e _v/ ^ 7r / 2 J_ iv ^/ 2 (pe J7r / 4 )), 0_ oc 0+ respectively. 
By adding linear combinations of <fi± of the form 

0+ + e 2i ^0_ (33) 

to the domain T>(t ) one obtains a one paramter family of self-adjoint extensions r 7 
of the operator r 0jOO . Solutions fl25|) of the eigenvalue equation belong to £>(t 7 ), if 
one chooses | = 2-^/ 2 r(l - iy/7/2) sin( 7 - ^f L )/r(l + i>/7/2) sin( 7 + ^). 
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